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FOREWORD 


Because  of  their  simplicity,  polytropic  models  have  historically  been  used 
to  provide  approximate  physical  representations  for  a  variety  of  classical 
astrophysical  objects  such  as  stars,  planets,  and  globular  clusters.  Since 
closed-form  solutions  exist  only  for  three  values  of  the  polytropic  index,  i.e., 
n  =  0,  1,  and  5,  much  attention  has  been  devoted  to  providing  analytic  represen¬ 
tations  that  are  valid  for  a  range  of  poly tropic  index.  This  work  was  conducted 
in  order  to  develop  convergent  power  series  solutions  for  the  radial  Emden  func¬ 
tions  ijjg  and  found  in  the  theory  of  rotating  polytropes.  The  practicality  of 
such  an  approach  is  that  when  used  in  astrophysical  applications,  mathematical 
operations  such  as  differentiation  and  integration  may  be  analytically  performed 
upon  these  functions  within  the  bounds  of  their  radii  of  convergence.  This 
report  was  reviewed  and  approved  by  Dr.  R.  J.  Anderle  and  Mr.  R.  W.  Hill. 

Released  by: 


THOMAS  A.  CLARE,  Head 
Strategic  Systems  Department 


ELECTE 
AUG  1  4 1985 


NSWC  TR  84-351 


CONTENTS 


Section 


Page 


INTRODUCTION. 


MATHEMATICAL  FORMULATION. 


POWER  SERIES  SOLUTIONS  FOR  ^Q(?)  AND  (5). 


CONVERGENCE  CHARACTERISTICS  OF  THE  AND  ^  SERIES 


CONCLUSIONS 


REFERENCES. 


DISTRIBUTION. 


•.*  *.*  */  %*  %V  *%*  vVV  VV  ‘  V/A^ •/  •.*  s*  -  M 


NSWC  TR  84-351 


1 .  INTRODUCTION 


In  recent  years  there  has  been  a  renewed  interest  in  the  study  of  polytropic 

gas  spheres  and  their  astrophysical  applications.  New  and  refined  numerical 

solutions  to  the  Lane-Emden  equation^  have  been  discussed  by  Seidov  and 
2  3 

Kuzakhmedov,  and  by  Suleiman.  Approximate  representation  formulae  for  these 

4  5 

solutions  have  also  been  developed  by  Service  and  by  Pascual.  Several  numerical 

studies  associated  with  the  characteristics  of  rotating  and  tidally  distorted 

poly tropes  have  been  reported  by  Naylor  and  Anand,^  Hachisu  and  Eriguchi,2 
8  9 

Jabbar,  and  Singh  and  Singh.  Convergent  power  series  solutions  to  the  Lane- 
Emden  equation  have  recently  been  developed  by  Seidov  and  Kuzakhmedov,^  and  by 
Mohan  and  Al-Bayaty . ^  * 

It  is  the  purpose  of  the  present  investigation  to  apply  the  method  of 
12 

Frobenius  to  the  rotating  polytrope  problem  in  order  to  obtain  power  series 
solutions  that  originate  at  the  center  of  a  polytropic  model  for  the  associated 
radial  Emden  functions  ^  and  In  section  2  the  equations  describing  the 

structure  of  a  slowly  rotating  polytrope  are  formulated.  The  power  series 
solutions  for  the  radial  Emden  functions  <|»g  and  are  developed  in  section  3, 
and  their  convergence  characteristics  are  discussed  in  section  4.  Conclusions 
drawn  concerning  the  present  study  are  presented  in  section  5. 


2.  MATHEMATICAL  FORMULATION 


Consider  the  equation  of  hydrostatic  equilibrium  for  a  rotating  spheroid: 


VP  =  pV(V  +  V') 


where  P  is  the  pressure,  p  the  density,  V  the  gravitational  potential,  and  V' 
the  rotational  potential  given  by 


V'  =  j  w2r2 [ 1  -  P2(m) ] 


1 
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Here  u>  is  the  angular  rotation  rate,  r  is  the  radial  distance  from  center, 

P2(u)  is  a  Legendre  polynomial,  and 

p  =  cos  z  (3) 

where  z  is  the  colatitude  in  spherical  polar  coordinates.  Taking  the  gradient 
of  Equation  (1)  gives 

^  VP]  =  V2(V  +  v')  (4) 

Using  the  Poisson  equation 

V2V  =  -4irGp  (5) 


where  G  is  the  gravitational  constant,  and  the  fact  that 

V2v'  =  2o)2  (6) 

in  Equation  (4)  gives 

V  [•—  VP]  =  -4trGp  +  2w2  (7) 

Assume  that  the  spheroid's  internal  structure  can  be  represented  by  a 
polytrope  of  index  n  so  that 

P  =  ICp1+l/n  (8) 

and 

j  p  =  Xan  (9) 

where  K  is  a  constant  related  to  the  mass  M  and  radius  R  of  the  configuration, 
and  X  is  its  central  density.  Since 

Vf-  ^P]  =  (n  +  1)K X1/rnV20  (10) 


(n  +  1)KX 


J.  /  u- 


4ttG 


1 


V2a  = 


n  _i_ 

-a  + 


Define 


2ttGX 


2 

0) 

2ttGX 


and 


C  =  $r. 

Then  Equation  (11)  becomes 

„2  n 

V  a  =  -a  +  a 


where  the  Laplacian  is  now  taken  with  respect  to  C  instead  of  r. 

Let  us  now  assume  that  a  solution  to  the  last  equation  to  first  order  in 
the  small  parameter  a  is  of  the  form 

a  =  0(C)  +  cx'F(C,p) 

where  0  is  a  solution  to  the  classical  Lane-Emden  equation 


dCU  dC ' 


+  0 


0 


(ID 


and  is  subject  to  the  boundary  conditions 
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After  substituting  the  relations 


2  2  2 
Vo  =  v^e  +  ctV  Y 


(19) 


and 


n 

°  ■  ~ e” 


+  an0n_1'l' 


(20) 


3=0 


into  Equation  (15)  and  equating  orders  of  the  parameter  a,  one  obtains 


V2f  =  -n0n  *¥  +  1. 


(21) 


Adopt  a  solution  to  the  last  equation  of  the  form 


*  *  +0(O  + 


00 


(22) 


k=l 


The  A^  coefficients  may  be  determined  from  an  examination  of  the  associated 
interior  and  exterior  gravitational  potentials.  Upon  integration  of  the  equation 
of  hydrostatic  equilibrium 


KX1/n(n  +  1 )Vo  =  V(V  +  V') 


(23) 


one  obtains 


.  .  ■  o  - 


int 


1  2 

-r  a£  [1  -  P»(p)I  +  constant 

D  Z 


(24) 


where 


int  K(n  +  l)X1/n 


(25) 


is  the  interior  potential.  The  external  potential  is  the  Laplace  solution: 
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j=l 


These  potentials  and  their  gradients  may  be  matched  at  the  surface 


6  -  i-e. 


Vt'V  ‘  W«l> 


$'  (£  )  =  $'  (£  ) 
int^l;  ext^l/ 


to  give 


1  _2, 


Q  (€ x)  +  a['l'0(51)  -  g  Cj]  ■  Y0/Ci  +  constant  j 
e'ap  +  -  3  Cj]  =  -Y0/Ci 


Al^ 1 1>  =  Yi /cj 


Al^l  (^1 )  =  _2yi^i 


A2*2^P  +  6  ?1  y2^1 


A2^2^P  +  3  ?1  3y2^1 


and  for  k  >  2 


1  k5 1 


-(k+1) 


Vi(ei>  =  -(k+  1)Ykrn 


- (k+2) 


(2 


(2 


(2 


(3 


(3 


The  conditions  of  Equations  (30)  define  Yq*  and  those  of  Equations  (29)  and  (31) 


imply  that  A^  =  0  (k  f  2).  is  obtained  from  Equations  (30),  which  yield 
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3.  POWER  SERIES  SOLUTIONS  FOR  ip  (£)  AND  i ^(O 


In  this  section,  the  method  of  Frobenius  is  applied  to  Equations  (38)  and 

(39)  in  order  to  obtain  power  series  solutions  for  the  radial  Emden  functions 

ip  and  i p  .  Development  of  the  series  from  the  center  will  only  be  considered 
u  1  10 

here.  Seidov  and  Kuzakhmedov  have  provided  such  a  development  for  the  Lane- 
Emden  function  0.  Since  their  results  are  used  in  this  work,  pertinent  expres 
sions  are  repeated  here: 


e-1  -  i 


where 


J 

=  ~  /  (in-j)a.c .  . 
J  i  J-i 


aj+l  *  j(]  +  l)‘(2j  +-3)Ijin  +  1  -  -  1  +  1)13 


+  2(j  -  i)la1aj+I.i 


c0  =1;  ax  =  -1/6 


and  j  >  1 . 


Consider  first  the  solution  for  Equation  (38) .  Let 


so  that 
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The  recursion  relation  for  successive  values  of  b .  may  be  obtained  by 
letting  j  -*  j+1  and  j  j-m  in  the  first  and  last  terms  of  Equation  (47), 
respectively,  giving 

°°  I  j-1 

7  I  2b .  , ( j  +  1) (2j  +  3)  +  nb.  +  n  7  b.  c 
JLmd  j  J+1  J  /  m 

j-1  (  m=l 

This  expression  has  been  equated  to  zero  since  there  are  no  vanishing  exponents 
of  ?  for  all  j  >  1.  It  should  also  be  noted  that  a  finite  upper  limit  has  been 
placed  upon  the  summation  index  m  since  j-m  >  1.  Thus  for  a  given  polytropic 
index  n,  one  has 
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2(J  +  l)(2j  +  3)  bj  +  7  bi-mCm 


The  cases  when  n=0  and  n=l  are  of  particular  interest.  Note  that  when 
n=0  one  sees  from  Equation  (50)  that  b^=0  for  j  >  1 ,  so  that 

2  12 
*0  =  V  =  6  ?  • 

When  n=l ,  0n  ^  =  1  and  the  summation  in  Equation  (50)  vanishes  giving 

,  .Ir2  _J_  ,4  1  f6 

’•'O  "  6  *  120  ^  5040  ^ 


,y>  (-i)J-1 

P0  ^_,<2j  +  1)! 


-2j  =  1  sin  £ 


These  results  thus  reduce  identically  to  the  known  closed  form  solutions  for 

g 

given  by  Jabbar  when  n=0  and  n*l. 

In  order  to  obtain  a  solution  to  Equation  (39),  let 


CO 

*2  -S, 


? -tj'V"" 


00 

■s 


2j(2j  -  l)d.^ 
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After  substitution  of  these  three  expressions  into  Equation  (39)  and  using 
Equation  (40) ,  one  has 


00 

00  ” 

B, 

J-l  1 

2 (2 j  +  3) (j  -  l)C2j~2  +  n 

1  +1}/" 

.  m=l 

523 

The  recursion  relation  for  d^  may  be  obtained  by  letting  j  -*  j+1  in  the 
first  term  and  j  j-m  in  the  last  term.  After  rearranging  terms  the  following 
expression  is  obtained: 


2j(2j  +  5)dj+1 


+  nd.  +  n  >  dJ  c 
3  /  J  j-m  m 

m«*l 


=  0 


( 


or 


-n 


J+1  2j(2j  +  5) 


+  >  d,  c 
/  j  j-m  m 

m=l 


for  a  given  polytropic  index  n.  Since  d^  is  arbitrary,  one  may  choose 


dl  =  L 


4.  CONVERGENCE  CHARACTERISTICS  OF  THE  ipQ  AND  SERIES 

The  convergence  properties  of  the  power  series  representations  for  and 
will  be  examined  in  this  section.  Specifically,  approximate  radii  of  conver¬ 
gence  will  be  established  for  these  series  by  comparing  results  obtained  from 
them  with  those  found  in  the  literature  that  were  computed  using  numerical 
integration  techniques.  Additional  comparisons  will  be  made  where  possible 
between  results  for  the  first  derivatives  of  the  series  and  their  numerically 
derived  counterparts. 
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Q 

Table  1  is  a  comparison  between  the  solutions  obtained  by  Jabbar  from 
numerical  integration  of  Equation  (38)  and  those  computed  from  the  series  given 
by  Equations  (44)  and  (50)  using  the  Naval  Surface  Weapons  Center  CDC  6700 
computer.  The  first  100  terms  were  used  in  the  series  evaluation,  and  the 
polytropic  index  n  was  allowed  to  range  in  increments  of  0.1  between  0  and  1.8. 
Both  sets  of  results  were  evaluated  at  £  =  £^,  where  the  £^  are  the  n  dependent 
first  zeros  of  the  Emden  functions  6.  Power  series  evaluations  have  not  been 
included  in  this  table  for  n  >  1.9,  because  the  series  diverge  for  this  range 
of  polytropic  index  when  £  =  £  ^ . 

As  can  be  seen  from  the  values  in  Table  1,  the  agreement  between  the 
numerically  integrated  results  and  those  obtained  from  the  100  term  series 
evaluations  is  quite  good,  with  the  level  of  agreement  generally  being  through 
three  to  four  decimal  places.  Increasing  the  number  of  terms  in  the  series  by 
50%  produces  negligible  improvement  for  all  n  considered,  whereas  reducing  the 
number  of  terms  by  50%  tends  to  reduce  the  level  of  agreement  with  the  numerically 
integrated  results  by  one  decimal  place  for  0.1  <  n  <  0.9. 

It  should  be  mentioned  that  the  polytropic  index  convergence  range  of  n  <  1.8 

for  the  4'n(£,)  series  is  consistent  with  the  convergence  range  of  n  <  1.5  when 

u  1  11  io 

£  <  £  1  determined  by  Mohan  and  Al-Bayaty  for  the  Seidov  and  Kuzakhmedov  series 

representation  for  0.  This  result  is  expected,  since  use  is  made  of  their  a^  and 
c^  coefficients  in  forming  the  b^  coefficients.  In  fact,  based  upon  recent  evalua¬ 
tions  of  the  0  series  made  by  the  author,  the  polytropic  convergence  range  for  0 
should  also  be  n  <  1.8  when  £  <  £^. 

Additional  comparisons  can  be  made  between  numerically  integrated  results 
obtained  by  Naylor  and  Anand^  and  Singh  and  Singh^  for  iJjq,  and  their  first 

derivatives  and  those  obtained  from  evaluation  of  the  first  100  terms  of 
Equations  (44),  (45),  (54),  and  (55).  These  results  are  presented  in  Table  2 
where  it  can  be  seen  that  the  agreement  is  quite  good  for  both  series  and  their 
derivatives.  The  level  of  agreement  is  somewhat  better  for  the  n  =  1.5  cases 
(especially  so  with  the  Naylor  and  Anand^  data)  than  for  the  n  =  2.0  cases.  This 
might  be  expected,  since  the  entire  interior  of  the  polytrope  is  contained  within 
the  radius  of  convergence  for  n  =  1.5,  whereas  the  values  of  ",  used  for  series 
evaluation  in  the  n  =  2.0  cases  are  nearly  the  same  as  the  radius  of  convergence. 

Estimates  of  the  radii  of  convergence  for  the  ip  and  ij>  series  were  made  by 
evaluating  these  series  at  increments  of  0.1£.  for  0  <  £  -•  :  and  increments  of 
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TABLE  1.  COMPARISON  OF  NUMERICALLY  INTEGRATED  AND 
SERIES  GENERATED  VALUES  FOR  *  (£  ) 


n 

?1 

(5  ^-Integrated* 

^0^) -Series** 

0.0 

2.44949 

1.00000 

1.00000 

0.1 

2.50454 

0.96500 

0.96586 

0.2 

2.56222 

0.93938 

0.94037 

0.3 

2.62268 

0.92205 

0.92288 

0.4 

2.68610 

0.91229 

0.91288 

0.5 

2.75270 

0.90963 

0.91002 

0.6 

2.82268 

0.91384 

0.91407 

0.7 

2.89628 

0.92484 

0.92496 

0.8 

2.97376 

0.94270 

0.94276 

0.9 

3.05543 

0.96764 

0.96768 

1.0 

3.14159 

1.00000 

1.00000 

1.1 

3.23261 

1.04024 

1.04021 

1.2 

3.32887 

1.08897 

1.08894 

1.3 

3.43081 

1.14693 

1.14690 

1.4 

3.53893 

1.21504 

1.21502 

1.5 

3.65375 

1.29440 

1.29438 

1.6 

3.77590 

1.38630 

1.38628 

1.7 

3.90606 

1.49231 

1.49230 

1.8 

4.04501 

1.61427 

1.61427 

results  pr 
t  100  term 


re  have  been  rounded  in  the  fifth  decimal  place 
n  used  to  evaluate  the  series 


ibar's 
le  firs 


esented  he 
is  have  bee 
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TABLE  2.  COMPARISON  OF  NUMERICALLY  INTEGRATED  AND  SERIES 


GENERATED  VALUES  FOR  *.,  *.,  AND  THEIR  FIRST 
DERIVATIVES  U  1 


Source 

u 

C 

*0 

3*0 

3C 

*2 

3*2 

3C~ 

1.5 

3.025 

Series^ 

0.95793 

0.46940 

4.12962 

1.08133 

S-S* 

0.95783 

0.46440 

4.12860 

1.08108 

3.20 

Series^ 

1.04232 

0.49659 

4.31342 

1.02551 

N-A** 

1.0423 

0.49660 

4.3134 

1.0255 

2.0 

3.47 

Series^ 

1.21098 

0.63453 

4.37464 

1.22236 

S-S* 

1.20729 

0.63147 

4.37357 

1.22206 

3.60 

_  .  t 

Series 

1.29604 

0.67434 

4.53587 

1.25675 

N-A** 

1.2961 

0.67495 

4.5360 

1.2617 

NSWC  TR  84-351 


5.  CONCLUSIONS 

Power  series  solutions  for  the  radial  Emdem  functions  and  have  been 
developed  at  the  polytropic  center  which  are  convergent  in  the  whole  interior  of 
a  rotating  polytrope  for  values  of  the  poly tropic  index  n  <  1.8.  Three  to  four 
decimal  places  of  accuracy  can  be  obtained  when  the  first  100  terms  are  used  to 
evaluate  these  series  in  this  region.  For  values  of  the  polytropic  index  n  >  1.8, 
the  radii  of  convergence  for  both  series  gradually  diminishes  until  n  *  3.3  where 
the  radius  of  convergence  is  approximately  0.1£p 

The  obvious  advantages  of  such  a  representation  are  that  within  the  radii  of 
convergence:  (i)  numerical  integration  of  Equations  (38)  and  (39)  is  no  longer 

required  to  evaluate  and  and  (ii)  mathematical  operations,  such  as 
differentiation  and  integration,  can  be  performed  analytically. 


15 
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